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A model to study two-proton emission from nuclei induced by electromagnetic probes is 
developed. The process is due to one-body electromagnetic operators, acting together with 
short-range correlations, and two-body A currents. The model includes all the diagrams 

■ containing a single correlation function. The sensitivity of the cross section to the details of 

the correlation function is studied by using realistic and schematic correlations. Results for 

\ the 16 nucleus are presented. 

-t— > 
O 

o 

m 



in 
o 
o 



- 1—^ 

X 



INTRODUCTION 



In this work we apply the model we have developed in these last years fH-|K to describe elec- 
tromagnetic responses by considering also short-range correlations (SRC) to the case of the elec- 
tromagnetically induced two-nucleon knockout from nuclei. The aim of our model is to take into 
account the SRC in electromagnetic processes involving nuclei with A > 4, in such a way that 
different processes, like inclusive electron scattering, one- and two-nucleon emission induced by 
electron scattering and, eventually, real photon scattering, can be described by using the same 
\ methodology. The idea was to have a consistent view of the different processes. 

. The linear response of the nucleus to an external operator O(q) can be written as: 



with 



S(q,u) = ImD(q,w) , 

7T 



£>(q,cj) = £ £n fa) (E n -E -uJ + ir])" 1 £„(q) 



where we have defined 



£n(q) 



(tt n |O(q)|tt ) 



This function involves the transition matrix element between the initial and the final states of the 
nucleus. These states are construct by acting with a correlation operator F on uncorrelated Slater 
determinants, as established by the Correlated Basis Function theory: 

|*o)=i ? |^o) |*„} = Fj$ n ). 
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Then the quantity £ can be written as 

($ n |F+O(q)F|$ ) 



fn(q) 



($ |F+F|<D ) 
We consider scalar correlation functions of the type: 

therefore the £ function results 

<*„|o(q) [[(i + MN 



($o|F+F|$q) 
($ n |F+F|$ n ) 



^n(q) 



($oin( i+ %)i*o> 



(* in( i+ Mi*o) 

i<3 

(*n|]I( 1 + MI*n> 

i<j 



where the the h function is denned as: 



hi 



1 



We perform the full cluster expansion of this expression and this allows us to eliminate the 
unlinked diagrams. At this point we insert the main approximation of our model by truncating 
resulting expansion such as only those diagrams which involve a single correlation function h are 
retained: 



&(q)-£(q) = <$n|0(q) i + 



i<j 



|$o)l 



In the above expression the subscript L indicates that only linked diagrams are included in the 
expansion. 

For one-particle one-hole final states, the function £ includes three terms: 



Cipih(q) = (*i P ih| O(q) l$o> + ($i P ih| O(q) E ^ l $ o) + <*i P ih| O(q) E h H l $ c 



The contributions of the various terms of the above expression are shown in Fig. ^ in terms 
Meyer-like diagrams. In addition to the uncorrelated transition represented by the one-point dia- 
gram (1.1). But, besides, also 4 two-point diagrams and 6 three-point diagrams are present. These 
new terms are necessary to get a proper normalization of the nuclear wave functions. In the fig- 
ure, the black squares represent the points where the external operator is acting, the dashed lines 
represent the correlation function h and the continuous oriented lines represent the single-particle 
wave functions. The letters h, i and k label holes, while p labels a particle. A sum over i and k is 
understood. 

Our model was tested by comparing our nuclear matter charge response functions with those 
obtained by considering the full cluster expansion 0. As we can see in Fig. |2]both calculations 
overlap and this gave us confidence to extend the model to other situations. 

With this model we have analyzed different processes involving both real and virtual photons. 
Inclusive (e,e') excitation of discrete nuclear states showed a very small effect due to SRC 1]. 
On the other hand, inclusive electron scattering in the quasi-elastic peak region presented SRC 
contributions smaller than final state interactions effects 2]. These last also dominate in the case 
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of (e,e'p) reactions 0]. Finally, in (7,p) we found a kinematic region showing pronounced sensitivity 
to the SRC. Unfortunately, in this region, the effect of meson exchange currents (MEC) is much 
larger than that of SRC [|. 

In all these processes the uncorrelated one-body response dominates the cross sections. To get 
rid of this contribution one has to investigate two nucleon emission processes. In this case one 
deals with two-body operators which can be either one-body operators acting in correlated nuclear 
states or MEC. If, in addition, one chooses two-proton emission, MEC involving charged mesons 
do not contribute and the situation is, at least in principle, rather clean to investigate SRC effects. 

In case one has two-particle two-hole final states, the function £ includes two terms: 

A A 

Ca P 2h(q) = ($2 P 2h|0(q) M $ o) + ($ 2P 2h|0(q) Yl M # o) 

Kj Ki<j 

The first term (see Fig. |3J) consists of 4 two-point diagrams. The second one results from the 
sum of 8 three-point diagrams. As in the previous case, this set of diagrams conserves the correct 
normalization. 

The geometry of the two-proton emission process produced by electron scattering from nuclei 
is shown in Fig. |3J The cross section is given by: 



dV ™ 2 



m pi P2 

" OMott free [VlWl + V t W t + V U Wtl + V tt W tt , 



de'dft e deidft pi dO p2 (2vr) 6 

where e' is the energy of the outgoing electron, Q e is the solid angle of its direction, ei is the energy 
of the one emitted proton, O pi and tt P2 are the solid angles of the directions of the two emitted 
protons, pi and P2 are their momenta, m is the proton mass, OMott is the Mott cross section and 
_f r ec describes the recoil of the residual nucleus. The t>'s are kinematic factors, while the w responses 
include the information about the nuclear structure and are given in terms of matrix elements of 
the charge and current operators. 

In our calculations we have considered the charge and the magnetization current one-body 
operators and the A-isobar two-body current, specifically these terms proportional to T3 0. 

The hole wave functions have been obtained from a Woods-Saxon potential whose parameters 
have been fixed to reproduce rms charge radii and single particle energies around the Fermi energy, 
while for the particles wave functions we used the optical potential of Schwandt et al. Q. 




FIG. 1: Meyer-like diagrams contributing to £ lplh (q)- 
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FIG. 2: Nuclear matter one-particle one-hole longitudinal response calculated with our model (full lines) 
and with FHNC (dashed lines) 
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FIG. 3: Meyer-like diagrams contributing to ^2 P 2h( c l) 



2. RESULTS 



The particular aspect we have investigated concerns to the sensitivity of the (e,e'2p) cross 
section to the details of the SRC. This has been done by using the three correlation functions 
shown in Fig. [3 The G (Gaussian) and S3 correlations, have been taken from a FHNC calculation 



© submitted to 6 th Workshop on "c-m induced Two-Hadron Emission" ,Pavia, 2003. February 8, 2008 - 5 




FIG. 4: Geometry of the (e,e'2p) process. 
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FIG. 5: Correlation functions used in our calculations. 



done with semi-realistic interaction . The V8 correlation is the scalar part of a state dependent 
correlation used in FHNC calculation done with a V8' Argonne interaction plus three-body Urbana 
IX interaction jjj]. 

The minimum of the G and S3 interaction is almost the same, while the V8 interaction has 
a deeper value. The V8 and S3 correlations overshoot the asymptotic value of 1 in the region 
between r=l and 2 fm. The cross sections obtained with these correlations for 16 O are shown 
in Fig. Apart the result of the 1 + state dominated by the two-body currents, all the other 
results have common trends. First, one should notice that the use of different correlations does 
not change the shape of the angular distribution. Second, the cross sections obtained with the 
Gaussian correlation, are larger than the other ones. 

To understand this result, we have done a set of calculations with rather schematic correlations. 
These correlations are shown in the lower right panel of Fig. The cross sections calculated with 
these correlations are shown in the other panels by lines of the same type. In these calculations 
the two-body A currents have not been included. 
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The box correlation indicated by the full line is our reference correlation. Lowering the min- 
imum (dashed lines) does not produce a large effect. The insertion of a part which overshoots 
the asymptotic value reduces the cross section (dotted lines). An analogous effect is obtained by 
reducing the size of the box (dashed-dotted lines). 

These results can be understood by remembering that the quantity entering in the cross section 
calculation is h(r) = 1 — f 2 (r) . The largest is the contribution of h to £^(q), the largest is the 
cross section. The overshooting of the asymptotic value, generates a term in h(r) of opposite sign 
with respect to the rest of the function, therefore the total contribution to the integral becomes 
smaller. The same effect can be obtained by reducing the size of the box as it is shown by the 
dashed dotted lines. 

To finish, we show in Fig. [8] the effect of each of the three ingredients of the transition operator: 
the two- and three-point diagrams coming from the one-body charge and current operators and 
the A-isobar current. As we can see, this last dominates the 1 + results, while for the 0^ the terms 
associated to SRC (in particular, these corresponding to 2-point diagrams) produce the larger 
contributions. In the case of the 2J, the 3-point diagrams are also relevant. 



We have presented here results for the (e,e'2p) process in 16 and 40 Ca obtained within a model 
which includes the effect of scalar SRC through diagrams with only one correlation line. 

We have found that qualitative features, such as the shape of the angular distributions, are not 
sensitive to the details of the SRC. On the other hand, slight differences in the correlation function 
produce large modifications in the size of the cross sections. 

However, the information about the SRC can be obtained only by a quantitative comparison 
between theoretical predictions and experimental data. Unfortunately the quantitative evaluation 
of the (e,e'2p) cross sections suffers from the uncertainties inherent to the required theoretical 
input. 

These kind of experiments are not the unique tool to extract the characteristics of the SRC 
correlations, but another procedure allowing, together with elastic, inclusive and one-nucleon emis- 
sion experiments, to obtain the information we are looking for. We are trying to describe all these 
process in a unique and coherent theoretical framework, thus permitting meaningful analysis of all 
of them. 
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FIG. 6: 16 0(e,e'2p) 14 C cross sections calculated with the correlation functions shown in Fig. The 
calculations have been performed for an electron incoming energy e = 800 MeV, a momentum transfer 
q = 400 MeV/c and a nucleus excitation energy ui = 100 MeV. The second proton was emitted with an 
energy e 2 = 40 MeV and within an angle 7 P2 = 60°. No recoil of the residual nucleus was considered and 
coplanar kinematics ((f> Pl = 4> P2 = 0) was assumed. The angular momenta refer to the state of the residual 
nucleus: (l P i/ 2 )- 2 : 0+ (lp 3/2 )- 2 : 0+ (lpi/2)- 1 {Ip^y 1 ■ 1+, (Ipi^rH^)- 1 ■ 2+ and (lp 3/2 )- 2 : 2+ 
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FIG. 7: Same as in Fig. El but for the schematic correlations shown in the lower right panel. 
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FIG. 8: 16 0(e,e'2p) 14 C and 40 Ca(e,e'2p) 38 Ar cross sections calculated with the same kinematics as in Fig. 
and with the Gaussian correlation. Dotted curves include only the 2-point diagrams of Fig. El Dashed curves 
include also the 3-point of the same figure (except the last four which so not contribute to this process). 
Solid curves show the results when also the A currents are included. 



